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Abstract:  Analytical  solutions  are  derived  for  a  one-dimensional  model  of  the  bulk 
temperature  response  of  open  channel  flow  with  unsteady  and  nonuniform  heating  at  an 
upstream  boundary,  the  water  surface,  and  the  riverbed.  The  model  describes  the 
temperature  variation  as  kinematic  waves,  and  the  solutions  are  explicit  formulas  that  are 
comprised  of  transient  terms,  which  play  dominant  roles  at  the  upstream  end,  and 
equilibrium  terms,  which  determine  the  temperature  far  downstream.  It  is  shown  that  the 
time-dependence  of  the  solutions  includes  solution  envelopes  that  can  be  both  complex  as 
well  as  interesting  because  of  the  interplay  between  the  upstream  and  lateral  boundary 
conditions.  The  applicability  of  the  solutions  to  practical  problems  is  demonstrated  for 
two  cases:  (1)  a  stream  bounded  at  its  upstream  end  by  a  dam  and  with  a  mid-reach 
inflow  that  is  also  subject  to  diurnal  heating;  and  (2)  Boulder  Creek,  Colorado,  which  is 
impacted  by  effluent  released  from  a  wastewater  treatment  plant.  The  model  prediction  is 
in  reasonable  agreement  with  gauged  data. 
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Introduction 


Although  accurate  and  computationally  affordable  numerical  simulations  are  becoming 
the  dominant  approach  for  predicting  mass  and  heat  transfer  in  various  flows,  analytical 
approaches  continue  to  be  used  because  of  their  irreplaceable  roles  in  many  studies  (e.g., 
Taylor  1921,  Fischer  et  al.  1979,  Tang  and  Sotiropoulos  1999,  van  Dongeren  and 
Svendsen  2000,  Chan  et  al.  2006).  This  paper  presents  a  model  for  the  bulk  temperature 
of  water  bodies  such  as  rivers,  lakes,  and  reservoirs  that  have  time-dependent  thermal 
forcing  at  their  upstream  ends  and  water  surface,  side  bank,  and  bed  (hereinafter  referred 
to  as  boundaries  ).  The  model  is  based  on  an  earlier  model  described  by  Edinger  et  al. 
(1974)  and  Jobson  and  Schoellhammer  (1987)  that,  in  addition  to  both  variants  and 
related  techniques,  has  been  applied  in  natural  streams  (e.g.,  Kim  and  Chapra  1997;  Boyd 
and  Kasper  2003;  Westhoff  et  al.  2007),  variable  reservoir  releases  (e.g.,  Carron  and 
Rajaram  2001),  and  river  temperature  control  (e.g.,  Gu  et  al.  1999).  A  comprehensive 
description,  including  supporting  parameters  for  phenomena  such  as  air-water  heat 
exchange,  can  be  found  in  Edinger  et  al.  (1974),  Jobson  and  Schoellhammer  (1987),  and 
Martin  and  McCutcheon  (1999).  It  should  be  noted,  however,  that  the  governing 
equations  of  the  model  are  typically  solved  numerically. 

This  study  considers  that  the  flow  is  fully  mixed  over  its  cross-section  and  that  diffusion 
is  negligible.  The  temperature  response  is  physically  simulated  as  a  series  of  one¬ 
dimensional  kinematic  waves  subject  to  unsteady  thermal  boundary  conditions.  The 
boundary  conditions  at  the  upstream  and  lateral  boundaries  are  either  periodic  or  arbitrary 
functions  of  time  and  space.  The  resulting  governing  equation  is  a  first-order  linear 
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differential  equation  that  can  be  analytically  solved.  It  consists  of  an  unsteady  term,  a 
convective  term,  and  a  source  term.  The  solutions  of  the  equation  are  direct  extensions  of 
the  equilibrium  solution  of  Edinger  et  al.  (1974),  which  is  valid  far  from  an  upstream 
boundary.  It  should  be  noted  that  analytical  solutions  for  similar  equations  have  been 
derived  in  other  studies  (e.g.,  Chatwin  1973;  van  Genuchten  and  Alves  1982;  Shukla 
2002;  Weigand  2004),  which  consider  longitudinal  diffusion  and  dispersion  of  a  flow 
with  an  initial  condition  and  unsteady  upstream  boundary  condition.  However,  these 
studies  did  not  consider  time-dependent  heating  and  cooling  along  the  channel.  Despite 
the  simplified  governing  equations,  the  solutions  presented  in  this  paper  can  provide 
estimates  for  practical  environmental  problems  such  as  the  construction  of  a  dam. 


Conceptual  Model  and  Governing  Equations 

The  model  is  depicted  in  Fig.  1,  which  illustrates  a  one-dimensional  flow  in  a  channel 
with  variable-temperature  water  entering  its  upstream  end  and  heat  input  along  its 
length.  Assuming  that  the  flow  is  well  mixed  over  its  cross-section  and  ignoring  diffusion, 
the  governing  equation  for  the  bulk  temperature  can  be  represented  using  the  principles 
of  thermal  balance  as  follows  (e.g.,  Edinger  et  al.  1974;  Jobson  and  Schoellhammer  1987) 


dT{t,s)  ,  .dT{t,s)  \  rrr  A 

— - - +  u{s) — - - =  K{T^^^{t,s)-T{t,s))  0<5<oo 

ot  os 

T{t,s)  =  T^^it)  s  =  0 


(1) 

(2) 
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where:  t  is  the  time;  5  is  the  distance  downstream;  T{t,  s)  is  the  temperature,  Tup  {t)  is  the 
temperature  at  the  upstream  boundary,  Tm  (t,  s)  is  the  lateral  boundary  temperature,  and 
u(s)  is  the  flow  velocity;  both  temperature  and  velocity  are  cross-section  averages.  The 
bulk  coefficient  of  heat  transfer,  .^^(defined  later),  is  dependent  on  wind  speed,  flow  depth, 
and  heat  capacity  of  the  water  body  (Edinger  et  al.  1974;  Jobson  and  Schoellhammer 
1987).  In  this  investigation,  is  a  constant.  The  time-dependent  temperature  at  the 
upstream  end  is  given  by  Eq.  (2),  and  the  unsteady  thermal  input  at  lateral  boundaries  is 
given  by  the  right  hand  side  (RHS)  of  Eq.  (1). 

Eq.  (1)  is  a  linear  partial  differential  equation  that  describes  the  propagation  of 
temperature  waves.  In  order  to  complete  the  general  problem  stated  by  Eqs.  (1)  and  (2), 
an  initial  condition  for  temperature  is  needed.  This  study  considers  only  equilibrium 
solutions,  which  are  independent  of  the  initial  temperature  condition, ,  and  thus  no  initial 
condition  is  used.  We  will  present  three  solutions  of  Eqs.  (1)  and  (2),  each  corresponding 
to  a  different  set  of  thermal  boundary  conditions  and  flow  velocity. 

Case  1:  Constant  velocity  flow  with  sinusoidal  time-dependent  thermal 
discharge  and  boundary  heating 

Case  1  has  the  following  flow  and  temperature  conditions 


u(s)  =  U 

(3) 

Tup(t)=%  +  ToSm((DQt  +  a) 

(4) 

=  T^+T^  sin  co^t 

(5) 
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where:  U  is  a  constant  velocity;  T^,  T^,  coo  ,  and  a  are  the  average,  amplitude,  frequency, 

and  phase  of  the  upstream  temperature  fluctuation,  respectively;  and  ,  T^,  and  coi  are 
the  average,  amplitude,  and  frequency  of  temperature  at  the  lateral  boundaries, 
respectively.  T^,  T^,  coo,  a  ,  T^,  T^,  and  coi  are  constants.  Conditions  (4)  and  (5)  are 

frequently  used  as  approximations  of  measurements  (e.g.,  Velz  1970;  Edinger  et  al.  1974). 
For  example,  Eq.  (4)  corresponds  to  the  temperature  modulation  of  water  released  from 
detention  ponds  and  reservoirs,  whereas  Eq.  (5)  represents  diurnal  heating  at  the  water 
surface.  Case  1  is  a  direct  extension  of  the  problem  discussed  by  Edinger  et  al.  (1974), 
who  used  Eqs.  (1),  (3),  and  (5)  to  describe  the  thermal  response  of  natural  water  bodies. 

Case  2:  Flow  with  spatially  and  time-dependent  boundary  heating 

Case  2  is  an  extension  of  Case  1  but  with  a  spatially  variable  lateral  thermal  boundary 
condition.  The  flow  velocity  and  upstream  temperature  are  represented  by  Eqs.  (3)  and 
(4),  but  a  spatial  and  time-dependent  temperature  is  used  at  lateral  boundaries 

'^bd  (^’  s)  =  T^+T^s  +  (Tj  +  Tj  s)  sin  co^t  (6) 

where:  and  are  the  downstream  gradients  of  average  temperature  and  amplitude, 

respectively.  Condition  (6)  can  represent  temporal  and  spatial  changes  of  atmospheric 
temperature,  riverbed  temperature,  groundwater  discharge,  and  snowmelt  at  higher 
elevations  (e.g.,  Hanrahan  2007,  Westhoff  2007). 
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Case  3:  Flow  with  spatially  variable  velocity  and  time-dependent  upstream 
temperature  and  boundary  heating 


As  another  extension  of  Case  1,  Case  3  incorporates  a  spatially  variable  flow  velocity  and 
time-dependent  temperature  at  upstream  and  lateral  boundaries 


u{s)  =  V(s) 

(7) 

II 

(8) 

Tbd  a)  =  71  +  S  +  p, ) 

(9) 

/=1 


The  right  hand  side  of  Eq.  (9)  defines  an  arbitrary  function  of  time  if  the  Fourier 
trigonometric  series  converges.  Therefore,  the  flow  velocity,  inflow  temperature  at  an 
upstream  boundary,  and  lateral  boundary  temperature  are  arbitrary  functions  of  distance 
or  time. 


Solutions  and  Discussion 


Case  1 


In  view  of  Lagrangian  coordinates,  T{t,s)  may  be  considered  as  the  temperature  of  a 
control  cross-section  by  letting  5  =  5(0 ,  5(0  being  the  traveling  distance  of  the  cross- 
section.  The  problem  represented  by  Eqs.  (1)  -  (5)  can  then  be  transformed  to 


dT(t,S(t)) 

dt 


KT{t,S{t))  =  KT,At,S{t))  0<f<oo 
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T(t,S(t))=T^^(t) 


(11) 


where:  to  is  a  reference  time,  and  dS(t)ldt  =  U  ,  or,  tQ  =  t  -  S(t)IU  .  Eqs.  (10)  and  (11) 
comprise  an  initial  value  problem  for  a  linear,  first-order,  non-homogeneous  ordinary 
differential  equation.  The  solution  can  be  obtained  as  follows  (see  Polyanin  and  Zaitsev 
1995) 


KZ 


T  (t,  s)  =  T^+  I  ^  =sin(fe)/  -  0^) 


To-Ti+T^sinico, 


t - 

V  Uj 


+  a>- 


KT, 


K^+o), 


rsm<  CO, 


t - 

V  Uj 


-0, 


•expl 


U  J  (12) 


where  0^  =  tan  '  (co^  /  K) .  Condition  (4)  has  been  used  to  eliminate  to  in  obtaining  Eq. 
(12).  Eq.  (12)  is  the  solution  for  Case  1  and  will  be  referred  to  as  the  base  solution 
hereinafter. 


The  base  solution  consists  of  two  parts:  (1)  the  first  two  terms  on  the  RHS  of  Eq.  (12), 
which  reflect  the  influence  of  the  lateral  boundary  condition,  are  the  equilibrium  solution 
far  downstream  as  discussed  by  Edinger  et  al.  (1974);  and  (2)  the  transient  solution 
represented  by  the  exponential  term  on  the  RHS,  which  results  from  both  the  upstream 
and  lateral  boundary  conditions.  The  two  parts  and  their  interplay  determines  the 
behavior  of  the  solution.  Only  the  first  part  dominates  the  solution  far  downstream 
because  the  second  part  is  damped. 
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The  Case  1  temperature  solution  given  by  Eq.  (12)  (Fig.  2)  oscillates  with  decreasing 
amplitude  until  the  equilibrium  region  is  reached  far  downstream.  The  solution  exhibits  a 
transient  stage  with  a  spatial  scale  that  is  about  the  same  whether  the  upstream 
temperature  is  higher  or  lower  than  the  equilibrium  .  Consequently,  a  characteristic 
transient  length  Lj  can  be  defined 


K 


(13) 


Lj  describes  the  effects  of  the  upstream  temperature;  the  transient  component  will  reduce 
to  5%  ( )  in  strength  at  a  distance  downstream  of  3Lt,  which  is  located  at  sKIU=3  in 
Fig.  2. 


The  solution  envelope  for  Case  1  (see  Fig.  2)  is  also  damped  downstream,  but  the 
damping  is  not  monotonic.  The  cause  of  this  behavior  can  be  explained  as  follows.  If  the 
frequency  of  the  upstream  temperature  oscillation  is  the  same  as  that  of  the  lateral 
boundary  temperature  ( =  co^)  the  base  solution,  Eq.  (12),  becomes 


T(t,s)  =  T^  +(Tq  -rj)-exp<{-^[>  +  Asin(fe)it-^i  +k) 


(14) 


where 
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where  A  is  the  amplitude  of  the  solution  envelope.  This  term  oscillates  downstream 
because  of  the  cosine  term  in  Eq.  (15),  but  it  is  slowly  damped  because  of  the  exponential 
term.  If  we  define  a  characteristic  envelope  length  as 

L,=—  (19) 
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it  can  be  seen  from  Eq.  (15)  that  the  envelope  wavelength  is  InL^ ,  which  is  determined 
by  the  velocity  and  temperature  frequency.  The  envelope  in  Fig.  2  has  a  spatial 
wavelength  of  lU  =0.84. 

An  interesting  situation  occurs  in  applications  such  as  river  temperature  control  when  the 
fluctuating  component  of  the  upstream  temperature,  as  given  in  Eq.  (4),  is  equal  to  the 

equilibrium  temperature  far  downstream;  Eg  =  KT^  /  ,  rUg  =  rUj ,  and  a  =  -0^ . 

For  example,  using  the  parameters  of  Fig.  2a  and  setting  the  upstream  temperature 
fluctuation  in  this  way,  the  resulting  solution  reveals  that  the  envelope  appears  as  parallel 
rather  than  wavy  lines  (Fig.  3).  In  fact,  the  base  solution,  Eq.  (12),  now  becomes 

T (t,  s)  =  T^+  I  ‘  sin(ruit  -  6>i )  +  [jg  -  Tj  •  exp 
^JK^+co^ 

The  amplitude  of  the  temporal  term  (second  term  on  the  RHS)  is  a  constant  and  thus  the 
extremes  of  the  solution  will  not  oscillate,  and  the  temperature  envelope  is  defined  by 
parallel  lines. 

Upstream  and  downstream  time  histories  (Fig.  4)  of  the  Case  1  solution  in  Fig.  2a 
demonstrate  the  superimposed  effects  of  the  upstream  and  lateral  boundary  conditions. 
The  resulting  upstream  temperature  oscillation  consists  of  multiple  frequencies  whereas 
the  downstream  solution  is  dominated  by  the  periodic  heating  and  has  a  single  frequency. 


|_&|  (20) 
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Case  2 


After  a  tedious  calculation  similar  to  that  for  Eq.  (12),  the  solution  of  Case  2  becomes 


T(t,s)=T,+T^s-  ^  sin(ruit  -6,) - f  ,  sin(rUit  -  6, ) 

K  K  +0), 


—  -  UT 

To-T,+  — 

0  1  ^ 


KUT,  . 
.2  .  , 


f  1 

r 

1 

1  kt,  .  \ 

[  1 

r 

t — 

+  a 

f  1 - ^ - sm( 

t — 

^  Uj 

J 

1  1 

^  Uj 

J 

(21) 


K^  +  a: 


t - 

V  Uj 


-e: 


•exp( 


& 

u 


1  2  2 

where  9^=  tan^  (2£Uj^  l{K  -co^  ) ).  This  solution  will  reduce  to  the  base  solution  of  Eq. 
(12)  under  conditions  (3),  (4),  and  (5).  The  solution  of  Eq.  (21)  with  the  flow  velocity  and 
upstream  boundary  condition  for  temperature  given  in  Eig.  2a,  but  with  temporally  and 
spatially  variable  lateral  boundary  heating,  is  plotted  in  Eig.  5.  Physically,  this  boundary 
condition  produces  an  increased  downstream  temperature  compared  to  Eig.  2a. 


Case  3 

rS(t) 

In  this  situation,  dS{t) !  dt  =  V {S{t)) ,  ox,  =  t  -  \  {I IV (r))dr  .  By  a  straightforward 

Jo 

derivation  similar  to  that  for  the  base  solution,  the  solution  is  obtained  as 
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(22) 

where  6^  =  tan  '  {co^  /  K) .  The  behavior  of  Eq.  (22)  is  illustrated  with  three  examples  that 
use  the  parameters  given  in  Fig.  2a,  but  with  different  flow  velocity,  upstream 
temperature,  and  lateral  boundary  condition.  The  first  example  (Fig.  6a)  demonstrates  the 
effect  of  a  spatially  variable  flow  velocity.  The  second  example  (Fig.  6b)  illustrates  a 
different  upstream  temperature  oscillation  and  Fig.  6c  employs  a  more  complex  lateral 
boundary  condition. 

Applications 

Open  channel  flow  with  time-dependent  discharge  and  surface  heating 

The  solution  of  Eq.  (21)  is  now  applied  to  a  stream  or  river  section  that  is  bounded  at  its 
upstream  end  by  a  dam  with  continuous  discharge  but  variable  water  temperature.  The 
stream  also  has  diurnal  heating.  The  river’s  mean  cross-sectional  temperature  can  be 
approximated  by  a  lateral  boundary  condition  that  is  not  only  a  sinusoidal  function  of 
time,  but  also  a  spatial  variable  because  of  a  change  in  air  temperature  along  its  length. 
The  temperature  at  the  dam  is  specified  as  having  the  same  frequency  as  that  of  the 
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surface  heating  but  with  a  phase  delay.  The  frequency  is  1/24  hr'^=  7.27x10'^  s'\  The 
river  just  downstream  from  the  dam  has  a  mean  flow  velocity  of  1  m  s'\  a  water  depth  of 
5  m,  a  mean  temperature  of  15  °C,  and  a  daily  fluctuation  of  3  °C.  The  air  has  an  average 
temperature  of  20  °C  at  the  dam  and  25  “C  at  500km  downstream  with  a  linear  increase, 
an  amplitude  of  15‘’C  and  a  phase  lag  of  -1.5  rad. 

This  problem  requires  calculation  of  the  heat  transfer  coefficient,  K,  which 
Edinger  et  al.  (1974)  derived  as 


K  = 


K 

pCph 


(23) 


3  — 

where:  p  is  water  density  (1000  kg  m'  );  K  is  the  air-water  exchange  coefficient  related 
to  factors  including  wind  speed  (100  W  m'^  °C'^);  Cp  is  the  heat  capacity  of  water  (4186  J 
kg-i  and  h  is  the  water  depth  (5  m).  Based  on  these  parameters,  K  equals  4.8x10'*^ 
s"'.  The  analytical  solution  is  plotted  in  Fig.  7,  with  the  upstream  and  lateral  boundary 
conditions  given  in  the  figure  caption.  As  estimated  by  Eq.  (19),  the  envelope  wavelength 
2  71  Le  is  86  km.  The  transient  length  3L7’is  obtained  from  Eq.  (13)  as  625  km,  suggesting 
that  a  river  with  these  characteristics  would  not  approximate  the  equilibrium  within  the 
distance  shown  in  the  figure. 


Consider  a  multiple  branch  flow;  suppose  the  channel  has  a  mid-reach  inflow  at  5  =  200 
km  that  raises  the  temperature  by  5  °C  and  reduces  the  average  flow  speed  to  0.5  m  s'\ 
The  solution  for  the  section  starting  from  the  location  of  the  mid-reach  discharge  to  its 
downstream  terminus  becomes  (s  >  200) 
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—  —  UT  T 

T(t,s)=Tj^+T^(s-  200) - -  +  I  ‘  =  sin(fe)i?  -  0^ ) 


K  Jk^+co" 


—  UT 
Hit)  +  5-T,+  ' 


KT 


f 


K 


^jK^+cof 


sin<  0), 


5-200 


t-- 
V  U 


-0, 


J 


(24) 

^(5-200)] 

r 


where  H{t)  equals  T{t,  200),  which  is  given  by  Eq.  (21).  Eq.  (24)  is  obtained  from  Eq. 
(21)  by  replacing  its  upstream  temperature  terms  with  [//(t)+5].  The  solution  of  Eq.  (24) 
is  also  shown  in  Eig.  7,  from  which  it  is  seen  that  the  solution  is  significantly  changed 
downstream  from  the  mid-reach  discharge  because  of  the  additional  inflow. 

Creek  flow  with  effluent  from  a  wastewater  treatment  plant 

The  analytical  solution  of  Eq.  (22)  is  applied  to  Boulder  Creek,  Colorado,  which  is 
impacted  by  effluent  discharge  from  a  wastewater  treatment  plant.  The  flow  is  shallow, 
highly  transient,  and  well  mixed.  A  complete  description  of  the  flow  and  a 
comprehensive  24-hour  survey  can  be  found  in  Windell  et  al.  (1988).  The  creek  is  13.7 
km  long  and  0.2  to  0.6  m  deep,  the  velocity  ranges  from  0. 12  to  0.4  m  s'\  and  the  creek 
temperature  varies  with  time.  The  air  temperature,  wind  speed,  solar  radiation,  and  other 
conditions  change  diurnally.  The  measurements  of  creek  temperature  were  made  at  four 
stations  located  at  0.6  km,  5.0  km,  9.0  km,  and  13.7  km  downstream  from  the  plant. 

Using  the  measured  flow  rates  and  creek  cross-sectional  areas  (Windell  et  ah,  1988),  it  is 
obtained  that  u  =  0.27,  3.10,  0.13,  and  0.13  m  s'*  at  station  0.6,  5.0,  9.0,  and  13.7  km, 
respectively.  Erom  the  survey,  it  is  also  determined  that  =  15  -r  10sin(7.27  x  10^^  t)  “C 
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{t  is  in  seconds),  K  =35  W  m'^  °C'\  and  h  =  0.4  m.  Here  is  obtained  by 

approximating  the  effects  of  air  temperature  and  solar  radiation.  Linear  distributions  of 
velocity  are  assumed  between  adjacent  stations  and  Eq.  (23)  is  used  for  the  heat  transfer 
coefficient.  Letting  station  0.6  km  to  be  the  upstream  end  and  approximating  the 

observed  temperature  with  T^p  =  17  +  3sin(7.27  x  lO^^t  -1.5)  “C  (Fig.  8a),  the  solution  of 

Eq.  (22)  for  the  mean  temperature  at  the  other  stations  is  given  in  Figs.  8b,  8c,  and  8d.  It 
is  seen  that  the  analytical  solution  is  a  reasonable  estimate  of  the  temperature  fluctuations 
measured  at  the  rest  of  stations;  the  analytical  solution  predicts  ranges  and  phases  for  the 
creek  temperature  similar  to  the  observations.  The  largest  difference  between  the 
prediction  and  the  measurement  occurs  at  furthest  downstream  station  (Fig.  8d). 

Although  the  temperature  predicted  by  the  analytical  solution  is  not  as  accurate  as  that 
from  a  numerical  model  by  Kim  and  Chapra  (1997),  the  errors  of  the  two  approaches  are 
of  the  same  magnitude  in  comparison  with  the  measurement.  The  results  are  promising 
given  that  the  flow  and  temperature  are  affected  by  many  factors  that  the  analytical 
solution  excludes,  such  as  the  material  comprising  the  creek  bed,  the  thermal  interaction 
between  water  and  sediment,  and  the  three-dimensional  character  of  the  flow. 


Concluding  Remarks 

This  paper  presents  a  one-dimensional  temperature  response  model  and  its  analytical 
solutions  for  flows  with  a  time-dependent  upstream  boundary  condition  for  temperature, 
and  temporally  and  spatially  variable  heating  along  lateral  boundaries.  The  solutions  are 
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expected  to  be  useful  in  providing  insight  into  the  physics  of  the  flows,  conducting 
parameter  studies,  and  setting  up  benchmark  solutions  for  testing  simulation  software.  As 
illustrated  in  the  previous  section,  the  solutions  produce  reasonable  estimates  of  the 
equilibrium  temperature  distribution  for  real  streams. 

The  analytical  solutions  are  valid  for  single  branches  of  rivers.  However,  as  shown  in  the 
first  example  of  the  previous  section,  they  can  be  applied  to  multiple  branch  rivers.  In 
order  to  account  for  flow  variability  and  boundary  conditions,  a  channel  can  be  divided 
into  several  sub-reaches  with  representative  values  for  the  temperature  variation  and 
heating;  the  solutions  can  then  be  applied  to  the  individual  sub-reaches.  The  analytical 
solutions  could  also  be  useful  in  estimating  contaminant  transport  problems  such  as  that 
in  a  stream  receiving  detention  water  (e.g.,  Wang  et  al.  2004).  In  addition  to  open  channel 
flows,  it  is  anticipated  that  the  solutions  are  applicable  to  heat  and  mass  transfer  in  other 
areas  (e.g.,  Morro  1977;  Maiani  2003),  as  long  as  they  are  one-dimensional  kinematic 
wave  phenomena.  The  analytical  solutions  are  explicit  formulas  consisting  of  simple 
functions  with  a  number  of  controlling  parameters,  but  the  calculations  of  the  solutions 
can  be  tedious.  We  propose  as  future  work  to  make  an  Internet  accessible  solution  tool, 
such  as  a  small  computer  code  or  Excel  spread  sheet,  which  allows  users  to  easily  set  up 
inputs  and  obtain  solutions. 

It  should  be  noted  that  the  analytical  solutions  are  designed  for  equilibrium  when  a 
balance  is  achieved  between  the  time-dependent  temperature  of  inflow  at  upstream  ends 
and  temporally  and  spatially  variable  heating  along  lateral  boundaries.  The  solutions  will 
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not  be  applicable  when  the  effects  of  the  initial  temperature  distributions  are  important. 
Moreover,  they  exclude  diffusion  effects,  which  are  often  important  in  practice.  In  this 
regard,  it  is  expected  that  the  solutions’  prediction  could  be  problematic.  For  instance, 
they  will  tend  to  predict  larger  temperature  extremes  as  diffusion  increases  in  importance. 
Also,  the  solution  accuracy  is  limited  because  factors  such  as  evaporation  and  the  three- 
dimensionality  of  the  problems  are  neglected.  To  include  these  factors,  a  numerical 
approach  must  be  adopted  (e.g.,  Jobson  and  Schoellhammer  1987;  Kim  and  Chapra  1997; 
Tang  et  al.  2008). 
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Notation 


The  following  symbols  are  used  in  this  paper: 

A  =  amplitude  in  Eq.  (14); 

B  =  parameter  in  Eq.  (15); 

Cp  =  heat  capacity; 

h  =  water  depth; 

Hit)  =  time-dependent  upstream  boundary  condition  for  temperature; 

Le  =  characteristic  envelope  wavelength; 

Lt  =  characteristic  transient  length; 

i  =  integers,  >  1; 

K  =  coefficient  of  heat  transfer; 

K  =  air/water  exchange  coefficient; 

5  =  distance  downstream; 

S( t)  =  traveling  distance  of  cross-section; 


t  =  time; 

Iq  =  reference  time; 

Tq  =  upstream  temperature  oscillation  amplitude; 

Tq  =  average  upstream  inflow  temperature; 

Tj  =  lateral  boundary  oscillation  amplitude  for  temperature; 

Tj  =  downstream  gradient  of  the  amplitude  of  the  lateral  boundary  condition 

temperature; 
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Tr 

T, 

T(t,s) 


Tbd(t,s) 


Tup{t) 


u{s) 

U 

V(s) 

a 


A 


(O^ 

fe)j 


0. 

o; 


Y,  K 


P 


=  amplitude  of  lateral  boundary  condition  temperature  due  to  component 

i; 

=  mean  lateral  boundary  condition  for  temperature; 

=  downstream  gradient  of  mean  lateral  boundary  temperature; 

=  mean  cross-sectional  temperature; 

=  lateral  boundary  condition  for  temperature; 

=  inflow  temperature  at  an  upstream  boundary; 

=  cross-sectional  mean  flow  velocity; 

=  a  constant  flow  velocity; 

=  a  variable  flow  velocity; 

=  phase  of  temperature  oscillation  at  an  upstream  boundary; 

=  phase  of  lateral  boundary  condition  oscillation  for  temperature  due  to 
component  i; 

=  frequency  of  temperature  oscillation  at  an  upstream  boundary; 

=  frequency  of  temperature  oscillation  at  a  lateral  boundary; 

=  frequency  of  lateral  boundary  temperature  oscillation  due  to  component 

i; 

=  phase  in  Eq.  (12); 

=  phases  in  Eq.  (22); 

=  phase  in  Eq.  (21); 

=  phases  in  Eqs.  (14)  and  (15); 

=  water  density. 
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Figure  Captions 


Fig.  1.  Schematic  representation  of  the  flow. 

Fig.  2.  Base  solution  of  Eq.  (12).  t/=  0.1  m  s'\  =50  °C,  and  K  =  0.04  s'\ 

r^^/^  =  l  +  0.4sin(0.30.  (a)r„/7;  =0.4  +  0.2sin(2t  +  1.5).  (b) 
r^^/7;  =1.6  +  0.2sin(2t  +  1.5) 

Fig.  3.  Base  solution  of  Eq.  (12)  with  the  upstream  boundary  condition  determined  by  the 
equilibrium  temperature  far  downstream.  t/=  0.1  m  s'\  =50  °C,  and  K  =  0.04  s'\ 

T^^IT,=  0.4  +  0.0326  sin(0.3t  - 1 .438) ,  /T)  =  1  +  0.4sin(0.30 

Eig.  4.  Time  histories  of  the  base  solution  of  Eq.  (12)  at  upstream  and  downstream 
locations.  t/=0.1  ms'\  Tj  =50  °C,  and =  0.04  s'\  =  0.4  +  0.2sin(2t +  1.5) , 

r^^/7;  =1  +  0.4  sin(0.30 

Eig.  5.  Solution  of  Eq.  (21).  U=0.l  m  s'',  7;  =50  ‘’C,  and  i^  =  0.04  s''.  /T)  = 

1  +  O.Obi’  +  (0.4  +  O.O65) sin(0.30  ,  and  7),^  /7j  =  0.4  +  0.2sin(2t  +  1.5) 
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Fig.  6.  Solution  of  Eq.  (22).  U=0.l  m/s,  =50  °C,  and  K  =  0.04  s'\  (a) 
u(s)IU  =  l  +  0.3sin(3.7Z5'/5) ,  0.4  +  0.2sin(2t +  1.5) ,  and  IT^  =  l  +  0.4sin(0.30  . 

(b)  u(s)/U  =1,  =  0.4 +  0.04(1 +  sin( 20)^,  and  =  l  +  0.4sin(0.30  •  (c)  u(s)/U  =1, 

T^pIT^  =  0.4  +  0.2sin(2t +  1.5) ,  and  =  l  +  0.4sin(0.30  +0.2sin(0.k -1) 

Fig.  7.  Temperature  response  of  open-channel  flow  with  a  mid-reach  inflow.  U  =\m  s'\ 
=15  +3sin(7.27xl0^'t-1.5)  °C,  =  20  +  10^^5 +  15sin(7.27 xlO^'O  °C, 

Tj '  =  10^^  °C  m'^  and  K  =  4.8x  10'^  s'\  t  and  s  are  in  seconds  and  meters,  respectively.  An 
inflow  is  located  at  5=200  km,  where  it  raises  temperature  of  the  flow  by  5  °C  and 
reduces  the  velocity  ioU  =  0.5  m  s-^(5>200  km) 


Fig.  8.  Temperature  response  of  a  creek  flow  with  effluent  from  a  wastewater  treatment 


plant.  Define  Error= 


-11/2 


/=1 


,  where  n  is  total  number  of 


measurements  during  the  24  hrs,  and  and  T‘  are  respectively  the  ith  measurement 

data  and  the  corresponding  predicted  temperature,  (a)  Station  located  at  0.6  km, 
Error=0.047.  This  location  is  set  up  as  the  upstream  end,  and  the  solid  line  in  the  figure  is 
actually  the  upstream  condition  of  the  analytical  model,  (b)  Station  located  at  5.0  km, 
Error=0.098.  (c)  Station  located  at  9.0  km,  Error=0.056.  (d)  Station  located  at  13  km, 
Error=0.13 
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